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1 Introduction 



In contrast to equilibrium systems, which are best analyzed using standard equi- 
librium statistical mechanics, there is no general approach to study the systems 
far from equilibrium. People are motivated to study the non-equilibrium sys- 
tems in one dimension, since these are in principle easier. Different methods 
have been used to study stochastic models in one dimension, including ana- 
lytical and asymptotic methods, mean-field methods, and large-scale numerical 
methods. Some models solved using these methods, are studied for example 
in SB- 

There is no universal meaning for the term exactly solvable. For example 
in |p^[T^|, solvability means that the evolution equation of n-point functions 
contains only n- or less- point functions. In solvability means that 

the S-matrix of the iV-particle system is factorized into products of 2-particle 
S'-matrices. This means that the .S'-matrices should satisfy the Yang-Baxter 
equation. Another meaning of integrability is that the time evolution equation 
for E n (t), the probability that n consecutive sites are empty at time t, is closed, 
that is it can be expressed in terms of other E m (t)'s. This method of solving 
the integrable models is called the empty interval method (EIM). 

The empty interval method has been used to analyze the one dimensional dy- 
namics of diffusion-limited coalescence |T^-|20|. Using this method, the functions 
E n (t) have been calculated. For the cases of finite reaction-rates, some approx- 
imate solutions have been obtained. EIM has been also generalized to study 
the kinetics of the (/-state one-dimensional Potts model in the zero-temperature 
limit 

In p2 ] , all one dimensional reaction-diffusion models with nearest- neighbor 
interactions, exactly-solvable through EIM, have been studied. In ]23], EIM has 
also been used to study a specific model with next-nearest-neighbor interaction. 
In |p4ft , the conventional EIM has been extended to a more generalized form. 
Using this extended version, a model not solvable by conventional EIM has been 
studied. 

In this article, we consider systems, in them more than two neighboring sites 
interact. We consider the most general systems with fc-site interactions. Some 
constraints are imposed on the interaction rates, so that the time evolution 
equation for E n (t) is closed. The general method of solving the time evolution 
equation is also discussed. Finally, as an example, a system with next-nearest- 
neighbor interactions has been considered in more detail. 



2 Models solvable through the empty interval 
method 

Consider a general one-species reaction-diffusion model on a one-dimensional 
periodic lattice with L + 1 sites, with a fc-neighboring-site interaction. We 
want to find criteria on the interaction rates, that guarantee the solvability of 
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the system via EIM, that is, the closedness of the evolution equation for the 
probability that n consecutive sites are empty, E n . 

Suppose that the initial condition of the system is translationally- invariant. 
Any configuration of k neighboring sites is denoted by a = (a%,a2,--- , ctfc), 
where ai — o or •. o(») is used to denote an empty (occupied) site. The rate 
of transition from a configuration a to b is denoted by A^. Similar to p2| , the 
interactions with k empty sites as initial or final configuration are not considered 
here. In other words for any a, 

A§ = A° = 0. (1) 

Excluding these interactions from the 2 k {2 k — 1) possible interactions, (2 fc — 
l)(2 fe — 2) interactions remain to be considered. We want to impose restrictions 
on A^'s in such a way that the evolution equation for l?„(f)'s be closed. As we 
will see, the form of evolution equation generally will be different for n > k — I 
and n < k — 1, and also will be different for n + k > L + 2 and n + k < L + 2. 
So we will treat each case separately. 



2.1 The case n > k - 1 and n + k < L + 2 

To obtain evolution equation for E n (t), one should first recognize the source and 
sink terms. There are two cases. In the first case, the intersection of the empty 
block and the interacting block is in the left-hand side of the empty block. In 
the other case, this intersection is in the right-hand side of the empty block. 
For the first case, the source terms come from 

n n 

a' i ■ ■ ■ a' i ci • • • c k l i o ■■■ o bi ^~o, (2) 

where c / 0. Here stands for a block of adjacent empty sites. One also has 
I < k — 1. A° = leads to I > 1. So the left source for E n is 

fc-l n-k+l 

^ = EEW c ^)- ( 3 ) 

1 = 1 a'.b 

Now consider the expansion 

a' ai SL2 a' 

=E£aLo, (4) 

/'=; a 

where in the last equality, a is an /'-dimensional vector and is (k — I' — 1)- 
dimensional. So, 

fe-lfe-l n+l-l'-l 

s * = E E E A -V(a • ^ )• (5) 

1=1 l'=l a ,b 
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If 

A«-=E A ^'-°.0' l<l<k-l, l<l'<k-l, (6) 

b 

n 

is independent of a, then one can sum up P(a • o • • • o) on the index a. Then 

fe-lfe-l n+l-l'-l 

1=1 i'=i 
fe-i fe-i 

= A "'^"+ i -''- 1 - E n+l-l>)- (7) 

i=l J'=J 

One can do similar calculations for the case that the intersection of the 
interaction block and the empty block is in the right-hand side of the empty 
block. Defining 

A«':=E A ^A,, l<l<k-l, l<l'<k-l, (8) 

b 

and assuming that it is independent of a, the source term for this case is 

fc-lfc-l n+l-l'-l 

^ = EE A "' p (^ •)• ( 9 ) 

i=i v=i 

Putting these together, the source term is 
fe-i fe-i 

S = E E( A »' + ^w)(E n+ i-v-i - E n+l - U ). (10) 
1=1 l'=l 

Now, lets consider the sink terms. Again we will treat interactions of left- 
and right-hand sides separately. First consider the left ones. The interactions 
which contribute to sink terms come from 

n n—k+l 

a'i • • • a' i o^~o bic x ■ ■ ■ c k _i 1T^~o, (11) 

where c ^ 0. The sink term from the left interactions is 

fe-i n 

i? L = -EE A aW^). (12) 
1=1 a',b 

Aq = leads to I > 1. One also has I < k — 1. is the transition rate, so 
it is defined only for a ^ b. But one can extend this definition and define the 
diagonal terms in such a way that 

E A ad=0- (13) 
b,c 



3 



Then it is seen that 

£ A aS = -E A aS- (14) 
b b 

Using this, one arrives at the following equation for Rl 

fc-l n 

^=EE A a'0^( a '^)- (15) 
1=1 a',b 

This again recasts to a simpler form, using the following expansion 

E] A a'i - a'jO = E/ A «'i -«i_i»0 E/ A ai - a;_ 2 »0 + ' ' ' + 
a' a' a' 

=EE Af aia ,.o. (i6) 

('=0 a 

Putting these together, one arrives at 

fc-l Z-l n+l-l'-l 

^ = EE A "' p (- ^ )> ( 1? ) 

where 

A«-=E A ".0' l<i<*-l, 0</'</-l, (18) 

b 

and it is assumed that Kf v is independent of a. 

It is seen that the conditions we have obtained for the source and sink terms 
for the left interactions, eqs. (^) and (|l8|), are similar, except for the range of 
I'. Performing similar calculations for the right interactions, all the conditions 
coming from the source and sink terms can be summarized as this. The following 
quantities should be independent of a. 

A^=Et«,!.o- 1<J<*!-1, 0<Z'<*-1, (19) 

b 

A w : =E A o.ai'\/> 1 <*<*-!, 0<l'<k-l. (20) 

b 

Defining A«/ := Afj, + Ajj,, the time evolution equation of E n (t), for n > k — 1 
and n + k < L + 2, takes the following form 

= EE A »'(^+'-"-i - E n+i-i>)- (2i) 

111 1=1 l'=0 

Note that in this equation, E is defined through 

E := 1. (22) 
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2.2 The case n < k - 1 and n + k < L + 2 

Now, we want to derive time evolution equation of E n (t)'s when n < k — 1. 
Two cases may occur. The first one, is that the n adjacent sites which we are 
focused on, are among the k interacting sites, and in the second case a block of 
these sites is outside those k sites. The result for the second case is similar to 
that of the preceding subsection, n > k — 1. For that case, we only quote the 
results. However, we study the first case in more detail. 
The source terms come from 



a\ ■ ■ ■ a„'ci . . . c„ei' • 



bi ■ ■ ■ bp o o ■ ■ ■ o di ■ ■ ■ d q , 



(23) 



where c ^ 0, p + q 



k, and p,q> 1. Then the source term is 
S= T £ A^ c d e ,P(a'ce'). 



P,3=l a',e',b,d 
p+q=k-n c ^o 



(24) 



Similar to the previous cases, one can rearrange the sum of the rates in the 
following form. 



a'.e' 



n— 1 q— 1 n— 1 

E E E A a«0 • • • O .e + E E A; 



f 

atOi 



(25) 



=0 n"=0 a,e 



=0 a 



where we have used the fact that c ^ and so at least one of the Cj's should be 



• . In the above equation, a is a (p 



l)-dimensional vector and e is a 



(q — n" — l)-dimensional vector. In the first term on the right-hand side of the 
above equation, the left • is the the first • in c from the right, and the right • 
is the first • in e' from the left. In the second term, the • is the the first • in 
c from the right, and e' is 0. Arranging all these together, one arrives at the 
following equation for the source term. 



E E 

p>9 =1 n'=0 



9-1 

E 

n"— a,e,b,d 



E E A a.°o.e^(a.0.e)+ ^ A^P(a.O) 



a.b.d 



Defining 



a _ V"^ \ 6i-"6pOdi-"d, 



• (26) 



(27) 



b.d 



and 



Pq,P' "~ / j A ai ■■■a r /»0 



b,d 



(28) 



where p<p'<p + n — 1, < q' < q — 1, p + q = k — n, and p,q>l. Assuming 
that A pq _ p / q > is independent of a and e and Ap g p , is independent of a, one can 
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sum up the terms in (Era 



s= E E 

p. <!=i n'=0 

or in terms of E n 's, 

n-l 

E E 

p,9=1 n '= 



9-1 



(29) 



9-1 



^pq,p'q' {E n i+ n " + -E"ri'+n"+2 — 2_E n /_|_,„//_|_l) 



(30) 



The independency of K pq ^ q i with respect to a and e, and / with respect to 
a is sufficient to guarantee that the above source term is expressible in terms of 
E n 's, but is not necessary. For example, in ( gij ) one can decompose the blocks 
c and a' instead of c and e', which leads to another set of sufficient conditions 
on the rates. 

Now, lets consider the sink terms for n < k — 1: 



'i ' ' ' "p " ' ' ' " c i ' ' ' e <j — * b\ ■ ■ ■ bpC\ ■ ■ ■ c n d\ ■ ■ ■ d q . 
The above interaction produces a sink term: 

R=~ E E A^P(a'Oe'). 



P. 3=1 a' ,e' ,b,d 
p+q=k-n c ^o 



Similar to the preceding case, using (0), one arrives at 
R = E E A a'°o d e^(a'0e') 



p>9 =1 a',e',b,d 

p+q=k~n 



Using the expansion 



9-1 p-i 

^a'Oe' ~ ^ai - a p /»0»ei - e g / 

a'.e' a.e g'=0p'=0 

9-1 P-I 

+ E E A 0. ei - e q , + E E ^ai -V* ' 

e q'=0 a p'=0 



(31) 
(32) 



(33) 



(34) 



G 



R can be written in the form 

*= E 



p+q—k— 



9-1 P-l 

kpq t pi q >(Ek-p'-q'-2 + Ek- 
<j'=Op'=0 



p'-q' 



2E k . 



p'-q'-l) 



9-1 

■ X A P9,9'(^ fe -9'-l ^ ^fc-g') 
?'=0 

p-l 

■ X/ ^pq.p'iEk-P'- 1 ~ E k-p') 
p'=0 



(35) 



where we have used the definition ( p7| ) and ( p8| ) for A- Pq y q ' and A£ p , but with 
an extension of the range of p' and q' to < p' < p + n — 1 and < q' < q — 1 . 
It has been also assumed that h. pq ^ p i q i is independent of a and e, and A^ , is 
independent of e. , is defined through 



P9.9' " 2-^1 0»ei - e t i 



pq-.q 

b± • ■■b p Od\---d q 



p + q = k — n, < q' < q — 1, 



(36) 



b,d 



and it is assumed that it is independent of e. Considering (|30|), (|35|), and the 
source- and sink-terms corresponding to the previous subsection, and noting 
that in these latter terms, one should replace 1 < I < k — 1 in the right-hand 
side of (|2l| ) with k — n < I < k — 1, one arrives (for n < k — 1 and n + k < L + 2) 
at 



dE n (t) 
dt 



k-l fe-1 



X E A-W (En+l-l'-l — E n+ i-i>) 



l=k-n l'=0 



p-l 



E E 

p-q =1 K. n'=0 

p+q=k—n 



9-1 

^ ] A pg,(p+ri-?i'-l)(g-n" — 1) 
n"=0 



X (£, 



n'-j-n" 



+ E n i +n " + 2 — 2E n ' +n 'i + i) 



A 



-n'-l(-^9+n' — -^g+n' + ly 



pq,p+n 



p-l <j-l 

^ ^ ^ ^ A-pq.n'n" (-^fe— n' — n ;/ —2 ~l~ Efc— n > — n " 2_£// i ;_ n /_ n '/_i) 
n'=0ra"=0 
p-l 

^ ^pq,n'( E k-n>-l - -Efe-n') 
n'=0 

E ^pq,n"i. E k-n»-l --Efe-n") 
ri"=0 J 



(37) 
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2.3 The case n + k > L + 2 

The last case to be considered is the case with n + k > L + 2. Normally the 
case of large L and finite k is of interest, in which one also has n > k. We 
assumed periodic boundary condition for the system. Then the intersection of 
the k interacting sites and the block of n sites may consist of two disconnected 
parts, of lengths I and So, one has, in addition to the source terms similar 
to those of subsection 2.1, a source term coming from 

a[ ■ ■ ■ a\ o • • • o b[ ■ ■ ■ b[, c\ • ■ ■ Odi • • • cZ fe -!-z' . (38) 



This leads to a source term 

n—l—l' 

S= E E A°?° a ,P(a'^b'c'). (39) 

i,!' = l a',b',c',d 
l + V = n + k- L-l .VO or b'^O 

Using 

z-i 



E ^i-^c'oi-oj ~EE A °_ 

a'b'c' p— C 

a'^O or b'^0 



• • o c»o • • • o 



i'-l 

+ • • • O »c o • • • o 

l-l i'-l 

+ EEE A °---°.c.o...o, (40) 

it is seen that if the quantities 

A ^, M :-E A °^-^2s ; 0<p<i-l, 0<g<i'-l 

d 9 V 

l' I 

A ^P : =E A |^c.|3, o< P <;-i 

i' i 

A; i ?, g :=E A ^-V^o, < g < Z' - 1 (41) 
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are independent of c, then the source term corresponding to (|38|) is 
( 1-1 1'-i 



S= E l E E Kl> ,pq(-2En+p+q-l-l' + l + E n+p+q -l-l> + E n+p+q -l-l, +2 ) 

l,l' = l ( p=0 q=0 

l+l'=n+k-L — l 

l-l 

+ ^2 Kl> ,p{En+p-l - -En+p-i+l) 
p=0 

I'-l } 

+ J2 Kv, q {E n+q -v - Wl'+O • ( 42 ) 

9=0 J 

Now let's consider the sink terms. Again there are terms similar to of sub- 
section 2.1, and a new sink term, which is 



r = - E E a'' : ; ..:/';a'^:. (43) 

l,L' = l a',b.c.d ^ Y 

l+l'=n+k-L-l b^O or d^O ! ' 

Using 

E ^ d = -E A a C0 (44) 



b,c,d 
b^O or d^O 



R can be written as 



R = E E A ° ::: oa°o-.-:°oP(a'^). (45) 

i,i' = l a'.c [5 ' J 

i + Z' = n + k-L-l 

Since a' 7^ (eq. ([I])), one has expansion 

k-i-i' 

E Ao • • • Oa' o • • • o = E E A ° ' ' ' o ai -a q »o- (46) 

a' """^ 9=0 a "^"^ 

If 



V I 



A ;^ : =E A ° ;:: °a. o'.'.'.°o, ;< P <fc-z'-i (47) 



is independent of a, then the above sink term becomes 

fc-i'-i 



R= E E ^U'A En +P- 1 ~ ^n+p-i+l)- (48) 



!,!'=1 p=i 
l+l'=n+k— L — 1 
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Note that here too, this condition on A' L is a sufficient condition for the EIM- 
solvability of the model. Using (f42|), (fl8|), and the source- and sink-terms cor- 
responding to those of subsection 2.1, one arrives at 



dE n 



L-n-l fc-1 



^ — ^ A;;/(E w+ ;_;/_x - E n+ i-i>) 

1=1 i'=0 

+ £ 



r f-i j'-i 



l,l' = l Lp=0 q=0 

l+l'=n+k—L — l 

fc-i'-i 

+ £ ^'iV ^ E n+p~l ~ -En+p+l-i) 
p=0 

i'-l 
9 =0 



(49) 



for n + fc > L + 2 (and n > k). Note that the summation limits in the terms cor- 
responding to the source and sink terms coming from the processes investigated 
in subsection 2.1, have been properly modified. 



3 General method of the solution 

In the previous section, the evolution equation of E n 's were obtained, eqs. (pl|), 
(|37|), and (p9|). Investigating ( |37j ) an d (|49| ), one can see that these equations can 
be rewritten in the general form of (|2 1[) , provided one defines E n 's for n < 0, 
and n > L + 1 properly. Doing this, one arrives at 



dE n (t) 
dt 



k-l fe-1 



^ ^ A;;'(£ ; „ + ;_;'_i — E n+ i-i>), 



(50) 



Z=l Z'=o 

for any n, with the following constraints (which are actually definitions). 

fc-i 

J2 M rs{E s - E s+1 ) = 0, -fc + 2<r<-l, 



(51) 



and 

r 

N rs (E s ^! - E s ) = 0, L + 2<r<L + fc-l. (52) 

s=L+2-fc 

In addition to these, there are two other boundary conditions 

E Q = 1, (53) 
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and 

E L+1 = 0. (54) 

This last condition comes from the fact that if the lattice is initially nonempty, 
it will never become empty (as it is seen from ([[])). So, excluding the empty 
lattice (which remains empty) there will always be at least one particle on the 
lattice. Equations ( |5l|) to ( ^4[) are 2k — 2 boundary condition for the difference 
equation (|50|), which is of the same order 2k — 2. To solve these equations, first 
consider the stationary solution. This solution (E^) satisfies 

fc-i fc-i 

E E ^U'{El +l -v-x J£h_,,) = 0, (55) 

1=1 l'=Q 

with the same boundary conditions ( pl| ) to (|54|). The solution to (|5^) is 

2fc-2 

El = E ( 5fi ) 

p=l 

where z p 's are the solutions of 

J2J2MAz l - l '- 1 -z l - l ') = 0. (57) 

1=1 l'=0 

This equation has 2k — 2 roots, one of them is 1. The coefficients a p can be 
determined using the constraints ( |5l| ) to (pi]): 

fc— 1 2fc-2 

E E M -«p( z p - z p +1 ) = 0. -k + 2<r< -1, 

s— r p— 1 
r 2fc-2 

E E ^.ctpC^- 1 - «') = 0, L + 2<r<£ + *-l, 

s=L+2-k p=l 

2fc-2 

E "p = 1, 
p=l 

2fc-2 

E ^ +1 - °- ( 58 ) 
p=i 

The full solution is of the form 

E n (t) =: E% + F n (t), (59) 
where F n (t) satisfies an equation similar to (pCj) but with homogeneous bound- 



11 



ary conditions: 



T7 = X] X/ hi' (Fn+l-V -1 - Fn+l-l'), 

m 1=1 l'=0 



fe-1 



53 M rs (F s — F s+ i) = 0, _fc + 2<r<-l, 

S — r 

r 

53 JV„(F a _i - F a ) = 0, L + 2<r < L + k-1, 

s=L+2-k 

F = 0, 

F L+1 = 0. (60) 
To solve this, one writes F n as 

F n (t) = J2e tt F e ,n, (61) 

e 

where F ei „ satisfies 

fc-i fc-i 



e^e,n = J! 53 ^■H'(-P , e,n+;-/'-l ~ ^e,n+J-J' ) , 
Z=l i'=0 

fe-1 

53M rs (F eiS -F e , s+1 ) = 0, -fc + 2<r<-l, 



53 A^i^s-i - F e , s ) = 0, L + 2<r<L + fc-l, 

s=L+2-fc 

= 0, 

F e , L+1 = 0. (62) 



F £i „ can be written as 



2fc-2 

Fe,n = E &, P < P , (63) 
p=l 



where z e ^s should satisfy 

fe-i fe-i 



5353W- r - 1 -^- r ) = e . (64) 



1=1 l'=0 
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This equation has 2k — 2 roots. The coefficients /3 e>p satisfy 

fc-1 2k-2 

E E M "&, P «P " <p') = °' -fc + 2 < r < -1, 

s— r p— 1 

r 2fc-2 

E £ Nrsfl^zl- 1 - < p ) = 0, L + 2<r<L + fc-l, 

s=L+2-fc p=l 

2fc-2 

E = °' 

P =i 

2fc-2 

p=i 

These are a set of 2k — 2 linear homogeneous equations for the 2k — 2 variables 
/3 e „. The condition that there exists a nonzero solution for these variables is 
that the determinant of matrix of coefficients be zero. This is a condition for e. 
So, in principle, one can solve this equation to obtain the solutions for e, and 
then the corresponding solution for z £l p's. One can then obtain /3 e , p 's, and F n {t) 
is obtained using ( |63| ) and (|6l|). 

4 A model with three-site interaction 

As an example, consider a model with three sites (next-nearest-neighbor) inter- 
action. Denoting the eight possible three-state configurations as following 

0:= (ooo) l:=(oo«) 2 := (o • o) 3 := (o • •) 
4:=(«oo) 5:=(«o.) 6 := (• • o) 7 := (• • •) , (66) 

and the transition-rate from the state i to the state j by , one can write the 
conditions for the solvability of the system through the empty-interval method 



13 



as 



4 


= Al, 














A 3 


= A|, 














A 4 7 


= Af = 


Af 


- A 4 










Ai 


= A^ = 


A^ 


= A 4 , 










A 2 


+ A 2 ! + 


A^ 


f A£ = 


\l- 


-Al 


+ Ai4 


" A6, 




+ A« + 


A§- 


f Af = 


A 2 7 ^ 


->§ 


+ A|M 


-A|, 


A? 


+ A? = 


A^ 


+ 4 = 


A^- 


f Ag 


= A 2 7 - 


f A?, 


A? 


+ A? = 


Ai 


+ 4 = 


Al- 


t- A| 


= At- 


f Ai 


A7 


+ A^ = 


A^ 


+ A 3 , 










4 


+ A* = 


A^ 


+ A 2 , 










A? 
















\ 2 

A 3 


- X 2 

— A ll 














\2 
A 6 


- X 2 















(67) 

For example, independence of Af^ with respect to a, gives A| = A|. One, of 
course, has also 

A* = A° = 0. (68) 

This is nothing but eq. ([!]). Using ([H]) forl = k- 2<n<L-k + 3 = L, we 
have 

1 = 1 l'=0 

= - A 20 E n+2 + (-A10 + A 20 - A 2 i)J5„+i + (A10 - An + A 2 i - K 22 )E n 
+ (An - A12 + A 22 )£n-i + A 12 -E„_ 2 , l<n<L. (69) 

The time-evolution equations for and £x come from (|37]) and (fl9|), respec- 
tively: 

d£?i(t) 



l'=0 



+ An,io(^o + ^2 - 2£Ji) + Af M (£i - E 2 ) 
+ A u , 00 (£i + £3 - 2£ 2 ) + A u>0 (£ 2 - S 3 ) 
+ Af li0 (^-B 3 ), (70) 



14 



and 

2 



dt 

l'=0 



+ A' llfi0 (E L ^ 2 +E L - 2B £ _i) + A£ i0 (£ L _! - E £ ) 

+ Ag^-Ei, - + A^oC^i-i - E L ). (71) 

These two equations can be rewritten in the general form of ( |69|) , provided one 
adds the boundary conditions corresponding to ([si]) and (|5^). These are in fact 
definitions of and £?l+ 2 : 
2 

5] A x ^(£i_^ - Ei-v) =A u ,io(£ + ^2 - 2^i) + Af Xil (£i - £ 2 ) 



i'=0 



+ A^o^ + £ 3 - 2£ 2 ) + Afi (£ 2 - £ 3 



+ Af 1)0 (£ 2 -£ 3 ), (72) 

and 



2 

A 2; -(i? £+1 _ i ' - E L+2 -i>) =A' 1100 (E L _2 + El - 2E L _i) 

l'=0 

+ A' 1 L lfi (E L ^-E L ) 

+ K L 1A (E L - E L+1 ) + Ai* (£ L _a - £ L ). 

(73) 

Equations (p9|), (|72"|), and (|73|) can be solved using the general method of the 
previous section. 

Now consider a special case 

A 7 = Ag = 0. (74) 
The conditions (|67|), and the nonnegativity of the rates, then lead to 
\\ = \\ = A? = 0, 

^2 = ^2 = 0' 

^3 = ^3 = ^3 — ^3 — ^3 — 0' 
-^4 = ^4 = ^4 = 0, 
^5 = ^5 = ^5 = ^5 — ^5 — 0' 
^6 = ^6 = ^6 = ^6 = 0' 

Ay - Ay = Ay — Ay = Ay = Ay = 0, (^5) 



and 



\5 _ \5 
A 2 — A 6l 

A 3 = ^2 + A 2 + A 2> 

\ 7 6 = \l + Xl (76) 
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Eq. ( p9| ) then reduces to 

E n = AE n+2 + BE n+1 -(A + B)E n , 1< n < L (77) 

where 

A:=Xl + Xl + Xl + \l, 

B := X\ + Xl + A| + Xl + Ag. (78) 

Eq. ( |70| ) becomes 

= A'Eg + 5'£ 2 - (A' +B')Ei, (79) 

where 

A' := Xl + 2Xl + 2Al + Xl + X 6 4 + 2X\ - X r 5 , 

B' := -Xl -Xl- 2Xl + A3 - Xl - Xl - 2X 7 4 + 2X\ + Ag, (80) 
and eq. ( |7l| ) becomes 

E L = B"{E L+1 -E L ), (81) 

where 

S" := Xl + Xl + Xl + X 3 2 + X 7 3 + X 5 4 + X 6 4 + Xl (82) 

This is in fact a degenerate example of the general case considered in the previous 
section. Note that E n = 0, 1 < £x+i is obviously a solution. This is expected, 
since the full lattice does not evolve, as A7 = 0. Noting that El+i = 0, one can 
solve ( |8l| ) to obtain El- This is found to be 

E L {t) = aL e- B,,t . (83) 

Using this, one can solve the equation for -Ex-i, to see that it contains two 
exponentials, exp(— B"t) and exp[— (A + B)t]. This is provided B" ^ A + B. 
(Note that in general B" < A + B. Equality holds iff A| = = 0.) Let us 
assume B" < A + B and proceed. It is not difficult to see that in other E n 's 
there are also terms like t l exp[— (A + B)t\. One can write 

L-n-l 

E n (t)=a n e~ B " t + /3nA~ (A+S) ', Kn<L + l (84) 

where 

a L+1 = /3 L+ltl = 0. (85) 
Putting this in j7^), one arrives at 

Aa n+2 + Ba n+l + (B" -A- B)a n = 0, (86) 



16 



and 



(/ + l)Au+i - A/3 n+2 j + Bf3 n+U . (87) 

The solution to ( pq ) is 

a n = a L ^ f , (88) 

Si - S2 

where ^'s are the roots of the equation 

(A + B - B")i 2 -B£-A = 0, (89) 



and ^ is arbitrary. The solution to (87) is 



1 m-i 



B A s 

P" A = £ (7— "tn+l+s, (90) 



s=0 



where 7 m 's are arbitrary constants if 1 < m < L, and zero otherwise. 

So far, all E^s except E\ have been obtained. Using (^), one can also 
obtain E\. It is seen that E\ contains similar terms and a new exponential 
term exp[— (A' + B')t\. So, in general there are only three time constants in 
the system, (as long as only the empty intervals are concerned). It may occur 
that two of these time constants, or all of them, are equal. This does not 
change the general behavior of the system. Only the degree of the polynomials 
multiplied in the exponentials are changed, and the corresponding coefficients 
can be calculated similarly. 
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